Increasing the Efficiency of a Thermionic Engine Using a Negative Electron Affinity 
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Most attention to improving vacuum thermionic energy conversion device (TEC) technology has 
been on improving electron emission with little attention to collector optimization. A model was 
developed to characterize the output characteristics of a TEC where the collector features negative 
electron affinity (NEA). According to the model, there are certain conditions for which the space 
charge limitation can be reduced or eliminated. The model is applied to devices comprised of mate- 
rials reported in the literature, and predictions of output power and efficiency are made, targeting 
the sub-lOOOK hot-side regime. By slightly lowering the collector barrier height, an output power 
of around lkW, at > 20% efficiency for a reasonably sized device (< lm 2 emission area) can be 
achieved. 



I. INTRODUCTION 

Producing more energy at high efficiency is a constant 
challenge facing humanity. Thermionic energy conver- 
sion has a role to play due to its many advantages: the 
design is simple, there are no moving mechanical parts, 
and there is no efficiency-reducing direct conduction of 
heat across the device 1,2 . 

Despite the advantages of thermionic engines, there 
have been technical challenges that have prevented adop- 
tion on a large scale. The challenges have been dealing 
with the high temperatures required to achieve a reason- 
ably good efficiency and having to overcome the effects 
of a space charged limited mode. 

The first challenge is to create a TEC which operates 
at a reasonably low temperature, but still outputs accept- 
able power at high efficiency In order to get the same 
electron emission from a material at lower temperature, 
one must either find a way to lower the thermionic barrier 
of that material, or one has to find a different material 
with a lower barrier. 

Naively, one might think that lowering the emission 
barrier of the emitter at all costs is the way to improve 
electron emission and therefore TEC output power den- 
sity. It is true that electron emission dramatically in- 
creases by lowering the barrier, but at the device level 
output power will not necessarily increase. One reason 
is that more emitter output current results in a greater 
negative space charge effect. Therefore emitter barrier 
lowering quickly hits a point of diminishing returns due 
to the self-limiting nature of the negative space charge 
effect. The second general issue stems from the fact that 
output power depends on the product of output current 
and the operating voltage. Negative space charge is- 
sues aside, the optimal operating voltage of the device 
is nominally the difference in barrier heights of the two 
electrodes. A strategy of lowering the emitter barrier 
without consideration of the collector will also quickly 
reach a limit of effectiveness as the value of the emitter 
barrier approaches that of the collector. 

The second challenge is the so-called negative space 



charge effect. The output current of a TEC is limited by 
the negative space charge effect because electrons travers- 
ing the intcrclectrode space create a negative charge bar- 
rier which blocks lower energy electrons from reaching 
the collector. Decreasing the interelctrode spacing is one 
approach to space charge mitigation. 

Generally speaking, the emitter barrier should be 
greater than the collector barrier. The emitter barrier 
should be low and its Richardson's constant should be 
high such that enough current is available to reach the 
desired output power. The electrodes should be close 
enough together to mitigate the negative space charge 
effect, bearing in mind that decreased spacing increases 
the engineering difficulty of device fabrication. 

Recent work has focused on improving the emis- 
sion characteristics of the emitter electrode. Nemanich, 
Koeck, and collaborators have reported favorable emis- 
sion parameters from nitrogen 3 and phosphorus doped 
diamond 4 . Further improvements in electron emission 
occur when a sample is irradiated with light via a phe- 
nomenon called photon-enhanced thermionic emission 
(PETE). Schwede et. al. first observed PETE from GaN 5 
and Sun et. al. have also observed it from diamond 6 . 

Regarding space charge mitigation, Hatsopoulous and 
Kaye demonstrated at considerable difficulty that a close- 
spaced configuration mitigates space charge 7 . More 
recently, Lee and collaborators have applied MEMS 
fabrication techniques to decrease the interelectrode 
distance 8 . Smith et. al. have predicted that an emitter 
which features a negative electron affinity, such as hydro- 
gen terminated diamond, could mitigate or eliminate the 
negative space charge effect 9 ' 10 . 

Improvement of the emission properties are welcome 
advances in the technology, but many materials already 
exist with quite good emission at low temperatures: scan- 
date electrodes, BaO, CsO, and other oxide materials 11 
all have exceptional electron emission properties. 

Little investigation is being conducted on improving 
the collector electrode technology. In this paper a model 
is developed which considers the effect of a negative elec- 
tron affinity collector on electron transport through a 
vacuum thermionic engine. It is shown that the nega- 
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FIG. 1. General schematic of a vacuum TEC. The emitter 
and collector electrodes are enclosed in a vacuum container 
and separated by some distance. The emitter is in thermal 
contact with a thermal reservoir at a higher temperature, and 
the collector is in thermal contact with a thermal reservoir at 
a lower temperature. Electrons are thermionically emitted 
from the emitter, travel across the interelectrode space, and 
are absorbed by the collector. The electrons travel from the 
collector through an electrical lead, through an external load 
where work is done, and back to the emitter to complete the 
circuit. 



tive electron affinity reduces, or in some cases eliminates 
the negative space charge effect. The model is applied to 
a TEC which is comprised of materials reported in the 
literature and could conceivably be constructed today. 
Finally, it is shown that of all the device parameters, the 
output power density and efficiency can be most effec- 
tively optimized by lowering the collector barrier. The 
model predicts that a reasonably sized device (< 1m 2 
emission area) can produce around lkW at > 20% effi- 
ciency at a hot-side temperature of w lOOOA". 



II. ELECTRON TRANSPORT 

A thermionic engine (aka thermionic energy conversion 
device or TEC) is a vacuum device that converts heat 
directly to electrical work. It is a heat engine and can be 
analyzed as such. A schematic diagram is shown in Fig. 
1. 

Electron transport is most easily understood by use 
of the electron motive and electron motive diagram; the 
motive diagram is similar to the band diagram in a solid- 
state device. Examples of a motive diagram in the case 




FIG. 2. Various conditions of electron motive in a TEC where 
the collector features NEA. The emitter is depicted on the left 
on each subplot, and the collector on the right. The maximum 
motive is indicated by a '+' symbol. 



where the emitter electrode has PEA and the collector 
has NEA are given in Fig. 2. Subscripts, E, and, C, de- 
note emitter and collector, respectively. The electrodes 
are separated by a distance, d. The motive at any point 
is denoted by -0, and the vacuum level at either surface 
is denoted by "0 appropriately subscripted (e.g. in the 
case of the emitter: ?Pe)- The maximum motive is de- 
noted by ip m . The quantity eV represents the potential 
of the collector relative to the emitter, V is referred to 
as the output voltage and e is the fundamental electronic 
charge. The output voltage is physically determined by 
the output current and the external electrical load of the 
device. For the case of PEA, <p denotes the thermionic 
barrier (sometimes referred to as work function), and in 
the case of NEA the barrier is denoted by ( which is the 
difference between the conduction band minimum and 
the Fermi energy. The value of NEA is denoted by \. 
The quantity tpc.CBM is the conduction band minimum 
of the collector, and (i denotes the Fermi energy. 

Electrons inside the emitter with energy greater than 
?Pe are thermionically emitted into the vacuum. The cur- 
rent density of the emitted electrons is given by Richard- 
son's equation. 



J = AT 2 exp - 



kT 



(1) 



where J is the emission current [Acm~ 2 ], A is Richard- 
son's constant, [Acm~ 2 K~ 2 ], T is the temperature [K], 
<f> is the barrier height or work function [eV], and k is 
Boltzmann's constant [8.62eV.K" _1 ]. 

Since the device is evacuated, electrons traveling across 
the interelectrode space represent a net negative charge. 
For some values of voltage, this negative space charge can 
cause the motive in the interelectrode space to be greater 
than the emitter barrier, presenting an additional space 
charge barrier to the thermionic electrons. This negative 
space charge effect decimates the output current for most 
devices and is a significant challenge in creating a viable 
thermionic device. 



III. THEORY 

Electron transport in the face of the negative space 
charge effect is modeled using a Vlaslov-Poisson system, 
following the success of this approach in the past 12,13 . 
Consider a TEC with a PEA emitter and a NEA collec- 
tor. We would expect this device should experience the 
negative space charge effect under most typical operat- 
ing conditions, but we assume the collector is sufficiently 
cool so that back emission is negligible. Electrons arriv- 
ing at the collector with energy greater than ipc,CBM are 
assumed to be absorbed by the collector. 

As the output voltage changes, this device will pass 
through several unique modes of electron transport. The 
accelerating mode is the condition where the maximum 
motive occurs at the emitter electrode. The saturation 
point is the greatest value of output voltage such that the 
maximum motive occurs just outside the emitter; it rep- 
resents a bound on the accelerating mode and is depicted 
in Fig. 2.1. 

As the output voltage increases, the maximum motive 
exists somewhere within the interelectrode space. The 
maximum motive barrier reduces the current traveling 
across the device, and so this set of voltages is known as 
the space charge limited mode. The motive diagram for 
a typical point in the space charge limited mode is de- 
picted in Fig. 2.2. The virtual critical point occurs when 
the maximum motive and collector conduction band min- 
imum coincide at the same height. This condition is 
shown in Fig. 2.3. 

As the voltage increases, the conduction band mini- 
mum of the collector limits the electron current entering 
the collector: electrons with energy less than the value 
of ipc,CBM cannot enter the collector and are scattered 
back to the emitter. 

At some value of voltage, the maximum motive occurs 
immediately outside the collector. This situation is re- 
ferred to as the critical point. For voltages greater than 
the critical point voltage, the TEC is in the retarding 
mode. 

Between the critical point and the virtual saturation 
point, the analysis of Langmuir 12 exactly models the elec- 
tron transport. For output voltage above the virtual sat- 
uration point, the electron transport differs from Lang- 
muir but the current is ultimately limited by the position 
of ipc,CBM- Therefore the output current is expressed by 
Eq. 2 which is unaffected by the precise details of the 
motive in the interelectrode space. 
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tion, me quantity, x m , is the position wnere ip m occurs. 
It is worth noting that the appearance of J in the expres- 
sion for £ is an indication that we are self-consistantly 
solving the Vlaslov and Poisson equations. These di- 
mensionlcss quantities allow Langmuir to write Poisson's 
equation in a universal form and thus numerically calcu- 
late the solution. Hatsopoulous and Gyftopoulous give 
a clear exposition of an algorithm to determine the out- 
put voltage given a value of space-charge limited current 
density 14 . Thus, the output current density vs. output 
voltage can be determined for a set of operating param- 
eters. 



A. Saturation Point 

In this case the output current is given by the satu- 
ration current of the emitter, and the saturation point 
dimensionlcss distance at the collector is given by 
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This value corresponds to a value of jes according to 
Langmuir's solution to the dimensionless Poisson's equa- 
tion. From Fig. 2.1 and the definition of 7, 



eV s = 



XC-Cc- lcskT E 



(4) 



B. Virtual Critical Point 



By the definition of the virtual critical point, and ref- 
erencing Fig. 2.3, 



XC = 1pm ~ ^CVR 



(5) 



The subscript CVR indicates Collector Virtual critical 
point. Using the definition of 7 and substituting, 



1CVR 



Xc 
kT R 



(6) 



Using the value of 7c va, one can use the dimensionless 
solution to Poisson's equation to determine £,cvr- From 
the definition of £, 



1 a^ I (c + eV 
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(2) 



Notation from Langmuir's analysis is used in the fol- 
lowing derivations and the highlights are listed here. 
In the space charge limited mode, Langmuir con- 
verts the motive and position to the dimensionless 
quantities; 7 = ^T 



is the dimensionlcss motive and 



I 



EVR 



6 



CVR 



-d 



/27rm e e^ 

V^p3" 



2\ 1/4 .1/2 



J, 



VR 



t: 



3/4 



(7) 



Using the dimensionless Langmuir solution, jevr can 
be determined from £,evr- From the virtual critical point 
motive diagram and the definition of 7, 



eVvR = <Pe - Cc + kT E jEVR 



(8) 



C. Space Charge Limited Mode 



The method of calculating the output current charac- 
teristic in the space charge limited mode is nearly identi- 
cal to the Langmuir case, one difference being the space 
charge mode extends from the saturation point only to 
the virtual critical point when the collector exhibits NEA. 
Hatsopoulous and Gyftopoulous's algorithm to find the 
output voltage given a value of output current density 
can be adapted to this case. 

1. Given J, calculate 7_e using Eq. 9. 

2. Calculate S,e using Langmuir's solution to the di- 
mcnsionlcss Poisson's equation. 

3. Compute t;c from Eq. 10. 

4. Again using Langmuir's solution to the dimension- 
less Poisson's equation, calculate jc from £q. 

5. The output voltage is given by Eq. 11. 
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IV. ANALYSIS 

A. The model applied to a plausible TEC 

Consider a TEC constructed from a scandate elec- 
trode described in 11 as the emitter, and a phosphorous 
doped diamond (PDD) electrode described in 4 as the 
collector. The emitter has a barrier of 1.16eV and a 
Richardson's constant of 7.8Acm~ 2 K~ 2 , while the col- 
lector has a barrier of 0.9eV and a Richardson's constant 
of 1 x 10 -5 ' AcvnT 2 'K~ 2 . No emissivity data is quoted 
for either material, but the scandate material is usually 
applied as a thin coating to a refractory metal such as 
tungsten. The diamond is CVD deposited as a thin film 
on molybdenum. A mid-range emissivity value of 0.5 is 
assumed for both electrodes. Efficiency of the device can 
be calculated by considering the electronic heat trans- 
port and the Stcfan-Boltzmann heat transport. No value 
is quoted for the NEA of the PDD collector. Accord- 
ing to the literature, the magnitude could be as large as 
leV 15 ' 16 , so a conservative value of 0.5eV is chosen. The 
interelectrode spacing is chosen to be 10/um. Finally, as- 
sume the emitter is held at a temperature of lOOOif while 
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FIG. 3. Output power characteristic showing the reduction 
to zero of the space charge limited mode by decreasing the in- 
terelectrode spacing. Inset: Motive diagram at the condition 
where the saturation and virtual critical point coincide. 



the collector is at 300K. At such low temperature, the 
back thermionic emission from the collector is negligible, 
consistent with the model. The output power character- 
istic of these parameters as calculated by the model is 
given by the dash-dot line in Fig. 3. The maximum out- 
put power density is 1.031Fcto -2 at an efficiency of 17.8% 
and occurring at an output voltage of 0.38V. Thus, the 
emission area of the device would be 0.1m 2 to produce 
lkW output power. 

The output power density/emission area of this device 
is reasonable, but the efficiency falls short of the 20% tar- 
get. Using the device just described as a starting point, 
individual input parameters can be adjusted (keeping all 
others fixed) to calculate the maximum efficiency. In 
this way, the parameters which have the greatest effect 
on device efficiency can be determined. Fig. 4 depicts 
the results of such a calculation as an array of plots. 
The first row of plots depicts the maximum efficiency vs. 
a particular input parameter; the columns from left to 
right show the effect of adjusting the emitter temper- 
ature, emitter barrier, collector barrier, collector NEA, 
and interelectrode spacing, respectively. The second row 
of plots depicts the output power corresponding to the 
maximum efficiency, and the third row of plots shows the 
voltage at which the maximum efficiency occurs. The red 
circles indicate the value of the common starting set of 
parameters. 

If there is a hard upper limit on the emitter temper- 
ature, the only way to achieve the 20% efficiency target 
is to decrease the collector barrier. At a value of 0.85eV, 
the collector barrier is sufficiently low that the device 
crosses the 20% efficiency threshold. This device pro- 
duces l.YJWcm output power density (corresponding 
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So Eq. 4 reduces to 




FIG. 4. First row: maximum efficiency vs. various TEC pa- 
rameters. Second row: corresponding output power density. 
Third row: corresponding output voltage. The red circles 
indicate the common starting set of parameters. 



to an emission area of 0.086m 2 for lkW total output) at 
an output voltage of 0.43V. 

The device can also achieve the target 20% efficiency by 
increasing the emitter temperature, ft is not surprising 
that increasing the temperature increases the efficiency 
because, in general, the carnot efficiency of a heat engine 
increases as the temperature differental increases. The 
emitter temperature must reach 1070^ before the device 
crosses the 20% efficiency threshold. 



B. Reduction/elimination of the space charge 
mode via collector NEA 



In some cases, the model of the TEC under considera- 
tion will experience no space charge limited mode. More 
precisely, there exists a surface in parameter space such 
that at the surface and points inside the surface, there is 
no range of voltages for which the TEC passes through a 
space charge limited mode. The initial device considered 
in this report will be on that surface in parameter space 
if the interelectrode distance is 3.89/im according to the 
following analysis. 

First, we derive the condition on that surface. The 
inset of Fig. 3 depicts the corresponding general motive 
diagram for a point on the surface. Consider a TEC with 
parameters on that surface: it is equivalent to say the 
saturation point and virtual critical point of that TEC 
coincide. We first compute the quantities at the satura- 
tion point. From the figure, 7cs becomes 
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(13) 



Note also that the output current density is equal to the 
emitter saturation current density. 

Next we compute quantities of the virtual critical 
point. From the figure it is clear that the output cur- 
rent density must also be equal to the emitter saturation 
current density. Moreover, £,evr equals zero and there- 
fore Jevr does as well. Eq. 8 becomes 
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The output voltage and output current density of both 
the saturation and virtual critical points are identical and 
therefore coincide. 

Another consequence of the fact that £evr = and 
Jvr — Jes is that Eq. 7 can be rewritten as 



£,cvr — d 
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Note that the left hand side of Eq. f 5 depends on the 
value of Jcvr via Langmuir's solution to the dimension- 
less Poisson's equation. If the above condition is met, 
the saturation point coincides with the critical point and 
the TEC is on the space charge limited mode boundary 
surface surface in parameter space. 

Notice also that there are device parameters on which 
the left hand side of Eq. f 5 depends, but for which the 
right hand side does not. For example, xc affects the 
value of £,cvr on the left hand side but no quantity on 
the right hand side. If a parameter is changed such that 
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(16) 



(12) 



then the device is operating within the space of parame- 
ters where it experiences no space charge limited mode. 
If the inequality is reversed, the device will pass through 
a space charge limited mode over some range of voltages. 



V. CONCLUSIONS 

A model was developed which considers the space 
charge limited electron transport through a TEC where 
the collector features NEA. Using this model, it is shown 
that NEA can mitigate, and in some cases eliminate, 
the negative space charge effect. Calculations were per- 
formed using the model to show output characteristics 
of devices made form available materials. Using material 
parameters quoted in the literature, and making reason- 
able inferences on missing parameter values, the output 



power and efficiency of a plausible device was calculated 
to be 1.03Wcm~ 2 and 17.8%, respectively. Since the ef- 
ficiency of the device was below the target of 20%, the 
input parameters were individually adjusted. The device 
can reach the 20% target when the collector barrier falls 
below 0.85eV. At that point, the output power density 
is l.nWcmT 2 . Such a device would output lkW with 
0.086m 2 emission area. Increasing the emitter tempera- 
ture to 1070-ftT would also achieve 20% efficiency. 



running OS X 10.8.2 (12C3103), Darwin kernel version 
12.2.1 with a 2GHz Intel core i7 and 8GB 1600 MHz 
DDR3 memory. The software to perform the calcula- 
tions in this paper was developed according to the best 
practices advocated by the Software Carpentry project 17 . 
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